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ABSTRACT 

The primary purpose of this thesis was to provide a 
method of checking the geometry and element connectivity 
input data for two finite element programs, ADINA and 
SAP IV. This preprocessor will accept the ADINA or SAP IV 
data deck, with minor modifications, and generate a graphi- 
cal display of the finite element model. The display is 
an oblique orthographic projection, and any orientation 
may be specified. Several options are available: exploded 

plots, partial plots, node numbering, element numbering, 
and others. Elements with three nodes on the same edge 
are plotted with a continuous curve on each edge generated 
by an interpolated parabola. Displacement postprocessing 
capability also exists. 
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I. 



INTRODUCTION 



A. GENERAL 

The continued developments and advancements of the 
finite element method this last decade have provided 
greater machine capabilities than ever before. Here at 
the Naval Postgraduate School (NPS) , two of the more 
favorable finite element programs are SAP IV [reference 1] 
for linear analysis and ADINA [reference 2] for nonlinear 
analysis. However, with the large amounts of numerical 
input/output data and automatic mesh generation, it is 
impractical to check and reduce this data without a graphi- 
cal representation. Data checking is divided into two 
areas: preprocessing and postprocessing. 

1 . Preprocessing 

Preprocessing is the checking of the input data 
deck. Errors in a finite element program occur basically 
in two areas. First, how close is the mathematical model 
(boundary conditions, loading conditions, material proper- 
ties, etc.) to the real problem? Second, are numerical 
errors present, or did misinterpretation of instructions 
occur in data deck preparation? Of the second type, most 
common errors are found in the geometry and element connec- 
tivity data. Preprocessing includes the forming of a 
graphical representation of the finite element model on 
which geometry and element errors are easily detected. When 
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node and element numbering options are available, the graph 
aids in the physical interpretation of the output data. 
Preprocessing is not a foolproof method of eliminating 
errors, but it does provide a tremendous advantage to the 
user. Preprocessors may be incorporated into the data check 
mode of the finite element program. However, to modify a 
large and complex program is dangerous. This may not be 
the best approach. A safer method is to develop a pre- 
processor which will read the finite element program deck 
separately with minimum modifications to that deck. 

2 . Postprocessing 

Though not as important as preprocessing, post- 
processing is extremely helpful in output analysis. Probably 
the most common and useful type of postprocessing is the 
contour plot. Appendix C of reference 4 lists a program 
developed to produce contour plots of stress data from finite 
element models. Contour plots can easily be adapted to a 
2D system, but 3D requires plotting the contours on 2D 
surfaces, a bit more complicated. Two methods of post- 
processing of displacements are the plotting of a deformed 
model or placing scaled vectors at the nodes. In the cases 
where the displacements are small, multiplication by a 
magnification factor produces an exaggerated representation. 
Like preprocessing, postprocessing can be incorporated into 
the finite element program directly, or done separately. 

When done separately, the finite element program must still 
be modified slightly to obtain a punched deck of the stresses 
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and displacements in the desired format. This thesis is 
primarily concerned with preprocessing. 

B. HISTORY OF DEVELOPMENT OF PSAPl 

A package [reference 4] containing digital computer 
programs for generating oblique orthographic projections 
and contour plots was produced by the National Aeronautics 
and Space Administration's (NASA) Langley Research Center 
(LRC) and distributed by the National Technical Information 
Service in January, 1975. The programs are completely 
general. Both programs contain options for selecting 
various plotting equipment including CALCOMP, VARIAN, and 
cathode ray tube (CRT) displays. With minor modifications, 
they can be adapted to most any system. 

1. SUBROUTINE PSAP Implementation 

Losh [reference 6], for his master's thesis in 
aeronautical engineering, implemented the preprocessor and 
postprocessor program, PSAP, at NPS in December, 1976. 
Modifying the LRC package [reference 4] , Losh adapted 
SUBROUTINE PSAP to the NPS IBM 360/67 system using the 
CALCOMP model 765 plotter. PSAP serves as a preprocessor 
for SAP IV models, and serves as a postprocessor for dis- 
placements of those models. Unfortunately, PSAP is severely 
limited in the type of elements it can plot. 

2. Motivation for SUBROUTINE PSAPl 

With the introduction of the ADINA [reference 2] 
program at NPS in January, 1977, and with expectation of 
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doing future analysis on ceramic turbine blades, it was 
desired to expand PSAP to include all ADINA elements and 
the 8-20 node brick elements in SAP IV. Like PSAP, PSAPl 
contains preprocessing and displacement postprocessing 
capabilities. PSAPl is presented in this thesis and has 
the following improvements over PSAP: 

a. Preprocessing for all ADINA elements. 

b. SAP IV 8 and 8-20 node elements. 

c. Expansion of SUBROUTINE ERROR. 

d. Interpolation of curves using shape functions 
[reference 3] through three points on the edges of the 
8-20 node brick elements and the 4-8 node plane elements. 

e. Improvements in defining the plot origin. 

f. Addition of an option (ISCALE =0) to plot 
sections of a model without losing perspective. 

g. Several other minor modifications considered 
improvements . 

C. PRESENT CAPABILITY 

PSAPl has the capability to plot all ADINA elements 
and all SAP IV elements except the pipe element. It will 
interpolate curves on the edges of 4-8 node plane elements 
and 8-20 node brick elements. Many options are given in 
Appendix A. Some of the more frequently used options are 
listed below: 

1. Numbering of grid points (NOTAT = 1) . 

2. Numbering of the elements (NOTAT = 2). 

3. Exploded plot (KDISP = 2). 
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4. Postprocessing of displacements (NUDISP or NVDISP 

or NWDISP = 1) in two forms: plot of deformed structure 

(KDISP = 1) or displacements represented by vectors at 
the nodes (KDISP = 3; see reference 6) . 

5. Symmetric representation about the XY (KSYMXY = 1) f 
XZ (KSYMXZ = 1) or YZ (KSYMYZ = 1) planes. 

6. Option to plot sections of the model (partial plot) 
to obtain a better view. Partial plots may be plotted to 
the scale of the complete model to avoid losing perspective 
(ISCALE =0) or blown up to obtain a better view (ISCALE = 1) . 
Multiple plots may be obtained using the same geometry and 
same displacement data (KODE = 1) , same geometry and new 
displacement data (KODE = 2), or new geometry and new 
displacement data (KODE = 3; see figure 1). 

In general, multiple plots (sections, partial plots, 
additional problems) present no problem. Plotting package 
user courtesy dictates that no more than 5 plots be placed 
on the CALCOMP plotter at any one time. Also, if the plots 
contain many elements (especially 8-20 node elements) , it 
is possible to run out of space in the plotting data sets. 

When this happens, you will receive 

ERROR IHC240I STAE, ABEND CODE IS: SYSTEM OB37 SYSPLOT. 

The best thing to do is split the run into two jobs. If 
the job must be run on one job (i.e., a large number of 
elements in the model or an assembly drawing where the 
scale of multiple plots is the same) , then SYSPLOT space 
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may be increased [references 7 and 8] by adding the card 
/ / GO . SYSPLOT DD SPACE= (CYL, (needed space) ), SYSOUT=C 
just prior to card 
//GO. FT10F001 DD UNIT = SYSDA 

in Appendix A. It would be wise to seek advice from a 
consultant in Ingersoll 146 if additional plotting space 
is required. 

D. EASE OF MODIFICATION 

Both PSAP and PSAP1 are written to maintain as much 
generality as possible for ease of expansion and modifica- 
tion. Several FORTRAN statements, variables and subroutines 
are not used. They were left purposely unchanged. Although 
PSAPl specifically reads ADINA and SAP IV data, it can 
easily be expanded to include any geometry and element data 
format. Simply study the read-in and storage methods (see 
Section II) , and construct appropriate subroutines to read 
any particular format. 
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II. PROGRAM ORGANIZATION AND DESCRIPTION OF OPERATION 



A. PSAP1 FLOW CHART 

Figure 1 is a condensed flow chart of PSAPl. Probably 
the most important information given on this chart is the 
sequence in which the data cards, NAMELIST OPTION and 
NAMELIST PICT are read. Remember, when generating a 
sequence of plots, once a parameter has been defined, it 
retains that value until it is reassigned. Note that when 
KODE = 1 or 2, the original values of NAMELIST OPTION and 
NAMELIST PICT are retained until they are changed. However, 
when KODE = 3 a new title card, NAMELIST OPTION and a set 
of problem data are read. All variables in NAMELIST OPTION 
and NAMELIST PICT are assigned their default values. A 
new problem begins in this case. It is important that the 
last NAMELIST PICT to be read must contain the value of 
KODE = 0. 

B. NAMELISTS AND EULER ANGLES 

1. NAMELIST OPTION 

Description and default values of NAMELIST OPTION 
are given in Appendix A. Basically NAMELIST OPTION varia- 
bles pertain to the given problem: the number of nodes, 

geometry format, displacement format, space between plots, 
and paper size. Default values for NAMELIST OPTION are 
set, and NAMELIST OPTION is read at the beginning of the 
problem. Assigned values will remain until exit from the 



16 




17 



program (KODE = 0) or a new set of problem data is read 
(KODE = 3) . 

2. NAMELIST PICT 

Like NAMELIST OPTION, NAMELIST PICT variable des- 
criptions and default values are given in Appendix A. 
Basically NAMELIST PICT variables pertain to a given plot. 
One NAMELIST OPTION may apply to several successive plots, 
but each NAMELIST PICT defines a unique plot. That plot 
may include the whole model, part of the model and any 
options defined in NAMELIST PICT. NAMELIST PICT also speci- 
fies the viewing plane through the Euler angles (figure 18, 
Appendix A) . 

a. Oblique Orthographic Projections (Euler Angles) 
An example of an oblique orthographic projection 
of a finite element model is given in figure 18 in Appendix 
A. The model can be viewed in any selected orientation. 
Euler angle transformations are used to specify orientation 
of the model to be projected. As described in reference 4, 
this transformation resolves the coordinate system of the 
model to a principal viewing plane (i.e., X q Y q , x 0 Z o' Y 0 Z o^ 
on which the display is to be plotted. Prior to rotation, 
the model coordinate system (X,Y,Z) is coincident with the 
coordinate system containing the viewing planes (X ,Y q ,Z o ). 
The viewing planes are fixed, and the model is rotated about 
its model coordinate system. The rotations (if»,0,< f>) of the 
body about the model axes (X,Y,Z) are shown in figure 18, 
Appendix A. The NAMELIST PICT variables KHORZ (horizontal 
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axis), KVERT (vertical axis), PSI (ip) , THETA (0), and 
PHI (<£) specify the viewing plane and Euler angles. The 
order of the Euler angle rotations must be PSI, THETA 



and then PHI. Mathematical transformations are: 
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SUBROUTINE ROTAT calculates the transformation matrices for 
every NAMELIST PICT, except if ISCALE = 0. ISCALE = 0 
directs the scale of the plot to be the same as that of the 
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previous plot. Should a rotation occur with ISCALE = 0, 
the plot width could exceed the paper width. 

b. Scaling 

The safest scaling method is automatic scaling 
(ISCALE = 1) . The user may specify a scale (ISCALE = 2) 
and the plot origin (XORGN, YORGN) , but one must be careful 
not to run the plotting pen off the graph paper. ISCALE = 0 
is a very useful option. The plot will use the same scale 
as the previous plot. It is useful in an assembly graph 
where examination of a mesh in sections without losing 
perspective is desired. Example 3, Section III, illustrates 
the option ISCALE = 0. When ISCALE = 1 in a NAMELIST PICT 
defining a partial plot, a "blow-up" of that section is 
obtained. ISCALE cannot be zero in the first NAMELIST PICT . 

c. Partial Plots 

To develop a partial plot, three methods of 
segregating elements exist: first, by the X, Y, and Z 

cutting planes; second, by node numbers, and, third, by 
element numbers. If a model has an area where the elements 
are relatively small, a "blow-up" may be desired. Choose a 
numbering scheme or coordinates to define the section to 
be segregated using one of the methods above. Example 3 
(figure 14, Section III) uses X, Y, and Z cutting planes to 
define the partial plots. Example 4 (figure 17, Section III) 
uses element numbers to section the plots. 
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C. NODAL POINT (GEOMETRY) INFORMATION READ- IN 

Nodal point data is read in by the GEOMn subroutines 
(GEOMl, GEOM2 , and GEOM9 , see figure 2). Since SAP IV 
and ADINA data decks are similar, SUBROUTINES GEOMl and 
GEOM9 are also very similar. They are both constructed to 
read and generate data in exactly the same way as ADINA and 
SAP IV. All data not needed by PSAPl is disregarded and 
the nodal point data is stored in array ZZZ (figure 3). 

After studying storage array ZZZ and GEOMl (or GE0M9) , 
a user moderately familiar with FORTRAN programming could 
easily construct a user supplied subroutine (GE0M2 ) to 
read the nodal point data in any desired format. 

D. ELEMENT (CONNECTIVITY) INFORMATION READ-IN 

After reading and storing the nodal point data, the 
element data must be read. The GEOMn subroutine will read 
the element control card (NPAR, references 1 and 2) . 

SUBROUTINE ELTYPE (figure 2) calls the proper element sub- 
routine to read the element data specified on the element 
control card. If several groups of elements are to be 
read, the process is repeated until all of the element 
groups have been read. Although the nodal point data is 
stored in array ZZZ, the element connectivity is read and 
stored on device 10 (disk) . 

E. DISPLACEMENT DATA READ-IN FOR DISPLACEMENT POSTPROCESSING 
Displacement data may be read in (figure 1) by SUBROUTINE 

DATA 9 ( KDATA = 9 ) or SUBROUTINES DATAl or DATA5 (user supplied. 



21 



! 1 




d 

-P 

d 

T 3 

>1 

-P 

•H 

> 

•H 

4J 

O 

CD 

G 

G 

o 

o 

4 J 

c 

CD 

s 

CD 

i— I 

CD 

-d 

g 

d 

>1 

u 

4J 

CD 

e 

o 

CD 

CT> 

cr> 

G 

d 

CD 

M 

J-l 

0 

M-J 

•P 

d 

,c 

o 

£ 

o 

I — I 

Pn 



CN 

CD 

D 

O' 



22 



ZZZ(l) ZZ Z (N+l ) Z ZZ ( 2N+1 ) ZZZ(3N+1) ZZZ(4N+1) ZZZ(5N+1) ZZZ(6N+1) 
ZZZ (2 ) Z Z Z (N+2 ) ZZZ (2N+2 ) ZZZ(3N+2) ZZZ(4N+2) ZZZ(5N+2) ZZZ(6N+2) 



z 

CM 

CSJ 

cm 





i — 1 


CN 




IS 
















1 " 


— 1 




"■ — • 














z 


Eh 


Eh = 


z 


Eh 














M3 


CU 


04 




04 
















s 


£ 




s 














CM 






















CM 






















CM 


i — 1 


CN 




a 
















"«— o' 


' — • 




" — • 
















Eh 


Eh = 


r 


Eh 
















CU 


04 




04 
















> 


> 




> 




1 








1 














«— 1 


CN 








m 












s 


S 






§ 


— ■ 


r— 1 


CM 


















CM 


> — • 


" — • 




> — ’ 














CM 


Eh 


Eh - 


z 


Eh 




1 — 1 


CN 






a 


cm 


a. 


04 




04 




> 


> 


** 




> 




D 


D 




D 


























1—1 


CN 








^ rtf 


t — 1 


CN 




a 




D 


D 






D 


Z — 


' — 1 


— 1 




— 1 
















Eh 


Eh r 


z 


Eh 


A 












— ' 


a. 


04 




04 


N — 


rH 


CM 






2 


CM 


CM 


Csl 




CM 




CM 


CM 






CM 


CM 






















CM 
























1 — 1 


CM 




a 




i — 1 


CN 






a 




— 


— 








X 


X 


* 




X 




Eh 


Eh = 


z 


Eh 














s 


04 


04 




04 














IS 


X 


X 




X 




1— 1 


CN 






a 


m 












X 


X 


** 


* 


X 


■ — ■ 




„ — ^ 




„ — ^ 














CM 


• — 1 


CN 




a 














CM 


— • 


— ' 










CN 


z 


z 


is 


CM 


Eh 


Eh = 


z 


Eh 




L_ 












04 


04 




04 
















X 


X 




X 




































IS 








„ s 














CM 


I— 1 


CN 




a 














v — • 


> — 1 






> — " 














CM 


Eh 


Eh 




Eh 














CM 


04 


04 




04 














Csl 


a 


a = 


z 


a 
















D 


D 




D 














^ — , 


1 !2 


a 




IS 1 


| 












IS 


L_ 





















CS3 

CM 

CM 



23 



Figure 3. Nodal point and displacement storage arrays. 

N = The number of nodes. (a) Array in subroutine PSAPl. 
(b) Arrays in subroutines called by PSAPl. (c) Nodal 
coordinates and displacements in (a) and (b) . 



KDATA = 1 or 5). When read, displacement data is stored 
in the last three columns of array ZZZ (figure 3) . PSAPl 
can postprocess displacements for both ADINA and SAP IV. 

The difficulty comes in obtaining a punched deck of cards. 
Reference 6 gives a description of how to obtain a deck of 
cards for SAP IV in a format acceptable to SUBROUTINE DATA9 . 
ADINA has no such provision. However, when preprocessing, 
the displacement data will be omitted (NUDISP=NVDISP=NWDISP=0) , 
and this step will be by-passed. 

F. PLOTTING LOGIC 

SUBROUTINE PLOTX (figure 4) is the main plotting routine. 
Since the nodal point data is stored in array ZZZ and the 
connectivity is stored on device 10, it is a simple matter 
to read the connectivity from device 10 (one element at a 
time) , and connect the nodes as they are defined in references 
1 and 2. For example, NEND = number of nodes defining the 
connectivity of a single element, NUMEL = the element number 
and NODE (NEND) is the array containing the connectivity. 

Device 10 contains this information successively for each 
element. The statement 

READ (10) NEND, NUMEL, (NODE (I) , 1=1 , NEND) 

will read the element connectivity to be plotted. The 
8-20 node brick (ADINA and SAP IV) and the 4-8 node plane 
elements (ADINA) may have 3 points defining each edge. If 
the midpoint node is defined , then isoparameteric 
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Figure 4. Flow chart for PSAPl plotting subroutine, PLOTX. 
Array NODE contains the connectivity of the 
element being plotted. NEND = the number of 
nodes defining the connectivity of a single 
element. NUMEL = the element number. 
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shape functions (reference 3) are used to interpolate along 
the three-node edges. These shape functions are identical 
to those used by ADINA and SAP IV, so the geometry represented 
graphically is identical to the problem solved in ADINA and 
SAP IV. SUBROUTINE CURVE does the interpolation with the 
following equations: 



X 

o 


= Nl * 


Xqi + N2 


* x 02 + N3 


* Y 

X 03 


Y 

O 


= Nl * 


Yqi + N2 


* Y 02 + N3 


* Y 

03 


Nl 


= S * 


(S - 1.0) 


/ 2.0 




N2 


= - (S 


+ 1.0) * 


(S - 1.0) 




N3 


= S * 


(S + 1.0) 


/ 2.0 





-1.0 <_ S <_ 1.0 

PSAPl uses the NPS plotting package [reference 5] 
subroutines . 
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III. PS API SAMPLE PROBLEMS 



The following examples have been chosen to illustrate 
some of the most useful options of PSAP1. Prior to 
attempting to use PSAPl, the user should have the problem 
defined and the cards punched in the format of references 
1 and 2. Appendix A of this thesis gives a complete des- 
cription of deck preparation for PSAPl here at NPS. This 
section should prove helpful in the understanding and 
interpretation of the options presented in Appendix A. 

A. ADINA EXAMPLES 

1. Reinforced Concrete Beam (example 1, figure 5) 

This example was chosen because it illustrates 
the importance of the exploded plot when more than one 
element group is used. It is taken from the ADINA manual 
[reference 1] . Figure 6 is a listing of the data cards 
as they are prepared for ADINA. Figure 7 indicates how 
that deck would be modified for PSAPl. Note: load cards 

are removed, NAMELIST OPTION and NAMELIST PICT are added, 
and the title to be plotted on the graph is added in figure 
7. Otherwise, figures 6 and 7 are the same. Figure 8, 
part (a) , illustrates an undistorted (KDISP = 0) PSAPl plot 
with the nodes numbered (NOTAT = 1). Figure 8, part (b) , 
shows the same mesh in an exploded form (KDISP = 2) with 
the elements numbered (NOTAT = 2). Note how the truss 
elements are visible in the exploded plot. 
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FINITE ELEMENT IDEALIZATION 



MATERIAL PROPERTIES : 

cr c = 3 740 psi 

<y \ - 4 58 psi 

°V Steel = 44000 psi 

E o = 6100 ksi 

concrete 

v = 0.2 

E 51^1 * 30000 ksi 
E . = 300 ksi 

T ste«l 



Figure 5. Example 1, ADINA 
truss and 8 node plane 
elements, Reinforced Concrete 
Beam (Given on page 84, 
reference 2) . 
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Figure 6. Example 1, ADIHA input deck listing, page 
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Not part of input deck. 
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Figure 7. Example 1, PSAPl input deck listing, page 
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Figure 7. Example 1, PSAPl input deck listing, page 
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2. Flat Plate With Hole (example 2, figure 9) 



This is a well known problem with which one can 
calculate the stress concentration on a hole in a plate 
under axial tension. Figure 10 is a listing of the PSAP1 
deck set-up. The mesh is composed of ADINA variable 4-8 
node plane elements. Notice on figure 11, parts (a) and 
(b) , how the interpolating shape functions round off the 
3-node edges. Part (a) has the nodes numbered (NOTAT = 1) . 
Part (b) has the elements numbered (NOTAT = 2) , and illus- 
trates the use of the symmetry option (KSYMXY=KSYMXZ=1) . 

The symmetry option enables one to obtain a picture of the 
complete plate even though the model only consisted of a 
quarter plate with proper boundary conditions. 

B. SAP IV EXAMPLES 

1. SAP IV Truss Problem (example 3, figure 12) 

Figure 13 is a listing of the PSAPl data deck. 
Figure 14 indicates how multiple partial plots can be used 
to obtain a better representation of the model. Part (a) 
of figure 14 is the complete model. Part (b) is the left 
half (XXMAX = 50' ), and part (c) is the right half 
(XXMIN = 50', XXMAX = 1.0E20'). Notice also that for 
Parts (b) and (c) , ISCALE = 0, which means succeeding plots 
have the same scale as the first. Had ISCALE equaled 1 
(blow-up) , then the width of the half view would have been 
the same as the complete model. Figure 14 size is limited 
by the NAMELIST PICT variable PLOTSZ. 
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Figure 9. Example 2 , flat plate with a hole in 

tension. ADINA 4-8 node plane element. 

Thickness = 1 inch 

Young's modulus = 30.0 X 10 6 psi 

Poisson's ratio = .333 

Total elements = 18 

Total nodes = 44 




Distributed Load = 

1000.0 psi 
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***** PSAP1 TITLE CARD TO OE PLOTTED ON GRA 

KIPLCR AE FLAT PLATF TESTING GEOM1 AND ADPLAN 18 ELEMENTS 
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Figure 10. Example 2, PSAPl input deck listing, page 
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Figure 10. Example 2, PSAP1 input deck listing, page 



***** NAMELIST PITT (ACTUAL STRUCTURE) 
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Figure 10. Example 2, PSAPl input deck listing, page 



Figure 11. Example 2, PSAPl output graphs, page 1 of 2. 




(a) Actual mesh, nodes numbered. 
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Figure 11. Continued, page 2 of 2. 
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Figure 12. Example 3, 
SAP IV truss element 
problem. 
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***** PSAP1 TITLE TO BE PLOTTED CN GRAPH 

SAP IV TRUSS EXAMPLE PROBLEM 
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Not part of input deck. 



Figure 14. 



Example 3, PSAPl output graphs. 




9 
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2. Cylindrical Bar With Spherical Hole (example 4, 
figure 15) 

This mesh could have several uses. Two might be 
to calculate stress concentrations if the bar is under 
axial load or to calculate loading if the void is under 
pressure (i.e., dispersed nuclear fuel pellet). Figure 16 
is a listing of the PSAPl data deck. Figure 17, part (a), 
is a representation of the complete model. Parts (b) , 

(c) , (d) , and (e) of figure 17 are partial plots of the 

total structure using options of the undeformed structure 
(KDISP = 0) with node numbering (NOTAT = 1) , and the 
exploded plot (KDISP = 2) with element numbering (NOTAT = 2) . 
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X fixed on YZ plane (X = 0) 

Y fixed on XZ plane (Y = 0) 

Z fixed on XY plane (Z = 0) 

Figure 15. Example 4, six inch 
a two inch diameter 
line under pressure 



Young's modulus = , 

30.0 X lOpsi 

Poisson's ratio = 

. 333 

diameter cylindrical bar with 
spherical void on the center 
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**** P SAP 1 TITLE TO BE PLOTTED CN GRAPH 
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Figure 16. Continued, page 4 of 6. ***** Not part of input deck. 
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Figure 16. Continued, page 5 of 6. ***** Not part of input deck. 



***** NAMELIST PICT FOR TCTAL STRUCTURE 
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***** Not part of input deck. 



Figure 17. Example 4, PSAPl output graphs, page 1 of 5. 
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Figure 17. Continued, page 2 of 5 . 




(b) Elements 1-10, nodes numbered. 



55 



Figure 17. Continued, page 3 of 5. 




(c) Elements 1-10, elements numbered, exploded plot. 
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Figure 17. Continued, page 4 of 5. 



113 115 113 123 127 




(d) Elements 11-20, nodes numbered. 
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Figure 17. Continued, page 5 of 5. 




(e) Elements 11-20, elements numbered, exploded plot. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 



With a little study of PSAPl and figure 2, one realizes 
how simple it would be to expand PSAPl to include virtually 
an unlimited number of geometry data input formats. Simply 
supply the GEOMn subroutines, element subroutines and 
modify SUBROUTINE ELTYPE if necessary. Another interesting 
project would be to incorporate PSAPl into ADINA and SAP IV 
so that a mesh plot could be obtained in the data check 
mode only. Reference 4 contains routines for plotting stress 
contours which could be incorporated. Hidden line logic 
is available. It would be difficult but could be incorporated. 
Hidden line plots are pretty, but they are not very practical 
when preprocessing because in preprocessing it is advantageous 
to see all of the nodes. 
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APPENDIX A 



PSAP1 USERS MANUAL 
I. NAMELISTS 

It may be useful for the user to review NAMELISTs in any 
good FORTRAN manual. However, a short description of the 
NAMELIST input is given here. Only columns 2 through 80 
of each card will be read. The computer expects to find a 
special delimiter symbol in column 2 of the first card 
followed by the NAMELIST name. The delimiter symbol for 
the IBM 360-370 series is the ampersand & ; other processors 
use the dollar sign $. Following the first card comes the 
list of variables and their assigned values. Some important 
points about the variables are: 

A. Variables and their values need not appear in any 
special order. 

B. Predefined variables need not appear in the list 
if it is not necessary to change the value. 

C. A comma should follow each assigned value. 

D. It is recommended that each variable appear on a 
separate card to facilitate change. 

The last card following the NAMELIST variables contains, 
beginning in column 2, the delimiter sign followed by the 
word END. For example, &0PTI0N (read NAMELIST OPTION) 
contains 12 variables. All are initially assigned default 
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values. But, suppose three of those values (NNDEST=200 
KGE0M=9, YSPACE=2.0) do not apply to a specific problem 
The &OPTION would be read as follows : 



&OPTION 
NNDEST=300 , 
KGE0M=1 , 
YSPACE=5. 0, 
SEND 



Blank space in column 1. 
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II. INPUT PROCEDURES 



The following sequence of cards is necessary to use 
PSAPl . 

NOTES IBM CARDS 
COLUMNS 

123456789 

(1) // [Standard green job card] 

// EXEC FORTCLGP, REGION. GO=150K 
//FORT. SYS IN DD * 

(2) C MAIN PROGRAM 

DIMENSION ZZZ(NZ) , DISPD ( 5 , 3 , NON) 

CALL PSAPl ( ZZZ , NZ , DISPD ,NON) 

STOP 

END 

(3) /* 

//LINK . USDD DD UNIT=3330 , VOL=SER=DISKO 2 , 

// DISP=SHR,DSN=S1153 .PSAPl 
//LINK. SYS IN DD * 

INCLUDE USDD (LOADM) 

ENTRY MAIN 
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( 4 ) 



/* 

// GO . FT10F001 DD UNIT=SYSDA, 

// SPACE= (CYL, (3,1) ) , 

// DCB= (RECFM=VS , BULKSIZE=3520 ) 



(5) // GO . SYSIN DD * 

(6) [PSAPl title card as it appears on plot] 

(7) &OPTION 

[&OPTION variables to be initially set or changed] 
SEND 

(8) [ADINA or SAP IV geometry data. Title 

card through element data — remove load cards] 

(9) [Case identification card — Omit if IDCASE = 0] 

[Displacement data cards — Omit for preprocessing 
only] 

(10) &PICT 

[ &PICT variables to be initially set or changed] 
SEND 

(11) [Additional data as defined by last value of 

KODE in NAMELIST PICT - Omit if last 
value of KODE = 0] 

( 12 ) /* 
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Notes : 



(1) Standard basic deck set up as described in chapter 
3 of reference 8. 



// [Standard green job card] 

// EXEC FORTCLGP, REGION. GO=150K 
/ /FORT . SYSIN DD * 



(2) Main program. 

DIMENSION ZZZ(NZ) ,DISPD (5 , 3 ,NON) 
CALL PSAPl ( ZZZ ,NZ ,DISPD ,NON) 

STOP 

END 



The main program has two functions: to allocate 

fast storage space, and to call PSAPl. NON must 
be greater than the number of nodes. NZ must be 
greater than 4*NON(7*NON if displacement data cards 
are to be input. i.e., NUDISP, NVDISP, or NWDISP = 1) . 



(3) Using load module library. 

/* 

//LINK. USDD DD UNIT=3330 , VOL=SER=DISK02 , 
// DISP=SHR, DSN=S1153 . PSAPl 
//LINK. SYSIN DD * 

INCLUDE USDD (LOADM) 

ENTRY MAIN 



PSAPl should be precompiled and stored in the machine 
in a load module. Load modules are a type of user 
library described in reference 7 and chapter 3 of 
reference 8. Since PSAPl requires over one minute 
to compile, precompiling results in a large time 
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savings. These cards are subject to change, and the 
most current version is determined by the user main- 
taining this library. If PSAP1 is not on a load 
module, these control cards may be replaced by 
subroutine PSAP1 (which consists of a box and a 
half of cards) . 

(4) Allocation of storage. 

/* 

//GO. FT10F001 DD UNIT=SYSDA, 

// SPACE= (CYL, (3,1) ) , 

// DCB= ( RECFM=VS , BULKS I ZE= 3520) 

PSAPl uses a slow storage device to store the element 
connectivity. These cards allocate 3 cylinders as 
described in chapter 3 of reference 8. 

(5) Deck set up card. 

//GO.SYSIN DD * 

This is a standard deck set up card described in 
chapter 3 of reference 8. 

(6) PSAPl title card. 

PSAPl title card is the title that will appear on the 
graph. Make sure a user identification is on this 
card. It consists of 80 alphanumeric characters. 

The first 40 characters will form the first title 
line. The last 40 will form the second line. 
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(7) NAMELIST OPTION. 



(Note: start in second column) 

&OPTION 

[&OPTION variables to be initially set or changed] 



SEND 




VARIABLE- DEFAULT 


DESCRIPTION 


NNDEST-200 


Must be equal to the number of nodes. 


NUDISP-0 


0 — X direction displacements not input. 

1 — X direction displacements input. 


NVDISP-0 


0 — Y direction displacements not input. 

1 — Y direction displacements input. 


NWDISP-0 


0 — Z direction displacements not input. 

1 — Z direction displacements input. 



(Note: unless displacement data is to be input. 



allow NUDISP, 


NVDISP , and NWDISP to default.) 


KGEOM-9 


Specifies the geometry input format. 

1- Subroutine GE0M1 reads in ADINA data 



deck geometry and connectivity. 

2- Subroutine GE0M2 may be supplied by 
user along with subroutines NSPLAN, NS3DEE and NSTRUS 
to read nodal data and connectivity in any format. 

9- Subroutine GE0M9 reads in SAP IV 
data deck geometry and connectivity. 

K DATA -9 Specifies the subroutine and corresponding 

method of input for displacement data. 



subroutine . 


1- Subroutine DATAl , a user supplied 



66 



2- Subroutine DATA 2 , a user supplied 



subroutine . 



9- Subroutine DATA9, reads a punched 
output displacement deck from execution of SAP IV as 
presented in reference 6. (Note: unless displacement 

data is to be plotted, allow to default.) 

NVALUS-0 Not incorporated, allow to default. 

IRESEQ-1 Not incorporated, allow to default. 

KPLOT-1 Not incorporated, allow to default. 

YSPACE-2.0 Space between plots in the Y direction 

in inches when successive plots are plotted (i.e., 

KODE t ^ 0). The graph title is plotted both on the 
top and at the bottom of each set of graphs controlled 
by a given NAMELIST OPTION. The space between the 
title and the plot is YSPACE/2.0. 

PSIZE-9.0 Paper size in the X direction, in 

inches. Used in scaling of the plots to insure this 
dimension is not exceeded. However, when manual 
scaling (ISCALE = 2; see NAMELIST PICT) this protection 
is not available, and it is possible to exceed the 
paper width. 



IDCASE-0 0- No identification card precedes 

the deck of displacement values. 

1- Identification card precedes the 
deck of displacement values . 
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(8) SAP IV or ADINA data cards. 



Here insert the geometry (node coordinates) and 
element connectivity. This includes the title card 
through the element data cards; the load cards are 
removed. For a SAP IV data deck, NAMELIST OPTION 
variable KGEOM = 9. For an ADINA data deck, 

KGEOM = 1 . 

From ADINA deck remove: 

1. Applied loads data. 

2. Frequency and mode shape calculations data. 

From SAP IV deck remove : 

1. Concentrated load mass data. 

2. Element load multipliers. 

3. Dynamic analysis cards. 

Otherwise, these cards are exactly the same as the 
deck prepared for SAP IV or ADINA. PSAPl is not 
limited to SAP IV and ADINA. The user may specify 
any unique format (i.e., KGEOM =2). In this case, 
SUBROUTINE GE0M2 must be supplied by the user. 

(9) Case ID card and displacement data cards. 

PSAPl is intended to be used essentially for prepro- 
cessing. If it is desired to use the postprocessing 
option, the user is referred to reference 6. Otherwise, 
omit the case ID card and displacement data. 



(10) and (11) NAMELIST PICT. 



(Note: start in second column) 



&PICT 

[&PICT variables to be initially set or changed] 
SEND 



VARIABLE- DEFAULT 



DESCRIPTION 



KHORZ-1 



Integer designating the horizontal axes 



of the viewing plane. 
(See figure 18 . ) 



1 = X . 
o 



2 = Y 



o 



3 = z 



KVERT-2 Integer designating the vertical axes 

of the viewing plane. 1 = X q . 2 = Y q . 3 = Z q . 
(See figure 18.) 



PHI-0.0 Angular rotation of the model about its 

X axis in degrees (performed 3rd, see figure 18). 
THETA- 0.0 Angular rotation of the model about its 

Y axis in degrees (performed 2nd, see figure 18) . 
PSI-0.0 Angular rotation of the model about its 

Z axis in degrees (performed 1st, see figure 18) . 
NEWFR-1 1- Frame change before plotting. 

2- No frame change before plotting. 
(Normally allow to default. A frame change resets 
the Y origin past the previous plot by YSPACE given 
in NAMELIST OPTION and resets the X origin at 0.0) 
ISCALE-1 0- No scale change. Use the same scale 

as the previous plot. Useful in an assembly graph 
where it is desired to examine a mesh in sections 
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page 127. 



without losing perspective. See example 3 in 
Section III of this thesis. ISCALE cannot be zero 
in the first NAMELIST PICT . 

1- Automatic scaling of plot and 
computation of proper origin location. 

2- User specified origin and scaling 

PLOTSZ-IO.O Maximum dimension desired on com- 

pleted plot prior to rotation. After rotation it 
is possible for the maximum dimension in the pro- 
jected plane to exceed PLOTSZ. The maximum size 
of the projection is unlimited on the vertical axis 
but is limited by the paper width (PSIZE in NAMELIST 
OPTION) on the horizontal axis. If rotation of the 
model causes the projection to go off the paper, 
it is rescaled prior to plotting. PLOTSZ is used 
in scaling only if ISCALE = 1) . 

XORGN-O.O X location of the plot origin. 

(Used only if ISCALE = 2.) 

YORGN-O . 0 Y location of the plot origin. 

(Used only if ISCALE = 2.) 

PSCALE-1. 0 Model size reduction factor. 

PSCALE is equal to the actual model size divided by 
the desired plot size. It is used only if ISCALE = 2 
(Note that when ISCALE = 2 is used, no rescaling 
occurs if rotation causes the projection to exceed 
the paper width.) 



NOTAT-O 



0- No numbering on plots. 

1- Numbering of grid points. 

2- Numbering of elements. 



XLHT-0.14 Height in inches of the integers 

specified by NOTAT. It should be a multiple of 
.07. If XLHT is not a multiple of .07, it will be 
rounded to the nearest multiple. XLHT has a 
maximum of .49 and a minimum of .07. 



KDISP-0 0- Plot of undeformed structure. 

1- Plot of deformed structure. 

2- Exploded plot. 

3- Displacement represented by 

vectors . 



KDISP = 1 or 3 represents a form 
of postprocessing and displacement data must be 
input in (9). If postprocessing is desired, refer 
to reference 6. 



IDMAG-2 1- Direct magnification of displace- 

ment data by DMAGS . 

2- Scaling of displacement data 
to a maximum value of DMAGS . 

DMAGS- 1 . 0 Magnification of displacements 

(if KDISP = 1 or 3) . Reduction factor of elements 
(if KDISP = 2) . 

KSYMXY-0 1- Symmetry about X-Y plane. 
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KSYMXZ-0 



1- Symmetry about X-Z plane. 



KSYMYZ-0 1- Symmetry about Y-Z plane. 

A plate quadrant with KSYMXZ and 
KSYMXY equal to 1 would yield a complete plate. 

See example 2 in Section III in this thesis. 

Note : To develop a partial plot, three methods of 

segregating elements exist: First, by X, Y and Z 

cutting planes; second, by node numbers; and third, 
by element numbers. The next ten variables are used 
to separate elements into partial plots. 



XXMAX, YYMAX, ZZMAX-1. 0E20 
XXMIN, YYMIN, ZZMIN- (-1.0E20) 



Local cutting planes 
Parallel to the principal 
planes . 



NDMAX-9999999 Maximum gridpoint identification 

number to be included in the plot. 



NDMIN-0 Minimum gridpoint identification 

number to be included in the plot. 



NELMAX-9999999 Maximum element number to be 

included in the plot. 

NELMIN-0 Minimum element number to be 

included in the plot. 
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KODE-O 



Specifies the control option after 



the plot is completed. 

0- Last plot, exit from program. 

1- Read another NAMELIST PICT. 

2- Read a new set of displacement 
data. (Postprocessing only.) For KODE = 2, dis- 
placement data must be followed by another NAMELIST 
PICT. 

3- Read a complete new set of input 
data starting with a title card. 

For KODE =1, 2 or 3 , additional 
sections of the deck must be repeated. The deck 
must end with a NAMELIST PICT having a value of 
KODE = 0. 

Note: A most important point to remember when 

generating a sequence of plots is that once a 
parameter has been defined, it retains that value 
.until it is reassigned. For example, if PLOTSZ is 
assigned a value of 8.0 for the first of a series 
of plots, and it is not redefined in any subsequent 
NAMELIST PICT; the value of PLOTSZ will be retained 
as originally specified. However, when KODE = 3 
and a new title, NAMELIST OPTION and problem data 
are read, all variables in NAMELIST OPTION and 
NAMELIST PICT are assigned their default values. The 
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problem starts over in this case. Refer to the 
flow chart in figure 1 and study the path for 
different values of KODE. 

(12) Delimiter card. 

/* 

Delimiter card is defined in chapter 3 of reference 

8 . 



III. SPECIAL FEATURES OF PSAP1 



A. POSTPROCESSING 

Reference 6 contains information and examples on the 
use of the postprocessor. SAP IV has the capability to 
punch displacement data cards in an acceptable format for 
PSAPl (specifically subroutine DATA 9 ) . As of this writing, 
ADINA does not have this capability, and data cards would 
have to be punched manually by the user. 

B. PARTIAL DATA 

Reference 6, on page 30, establishes a procedure by 
which it is possible to input only a portion of the finite 
element model for a data check. This feature is valuable 
in a case where several different people are preparing 
different parts of a large data base for a problem and 
desire to check individually their inputs graphically for 
accuracy. All nodal coordinates for the entire model may 
be input, or only those that specifically define the portion 
of the finite element model to be plotted. In either case, 
all nodal coordinates that define the elements to be plotted 
must be specified . To use this feature the "element control 
cards" (described in references 1 and 2) must be modified. 
For example, if only the connectivity for elements 15 
through 50 of a problem are available, the changes below 
would be made. 
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1. All SAP IV Elements 

All SAP IV element control card changes would 
have similar changes because columns 6-10 contain the 
total number of group elements. Columns 66-70 are not 
used. Thus to plot only elements 15 through 50, make the 
following changes: 

a. Enter the upper bound (i.e., 50) in columns 

6 - 10 . 

b. Enter the lower bound (i.e., 15) in columns 

66-70. 

2 . ADINA Truss, 2D and 3D Elements 

ADINA element control cards for the truss, 2D 
continuum and the 3D continuum elements would be changed 
as follows: 

a. Enter the upper bound (i.e., 50) in columns 

5-8. 

b. Enter the lower bound (i.e., 15) in columns 

53-56. 

3 . ADINA Beam Element 

ADINA element control cards for the beam element 
would be changed as follows: 

a. Enter the upper bound (i.e., 50) in columns 

5-8. 

b. Enter the lower bound (i.e., 15) in columns 

65-68. 
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APPENDIX B 



SUBROUTINE PSAP1 LISTING 
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